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Abstract 


Spherical harmonic analysis was used to analyze the observed clima- 
tological (C) fields of temperature at 850 mb, geopotential height at 
500 mb, and sea level pressure. .The spherical harmonic method was also 
applied to the corresponding "model climatological" fields (M) generated 
by a general circulation model, the "GISS climate model." The climate 
model was initialized with observed data for the first of December 1976 
at 00. GMT and allowed to generate five years of meteorological history. 
Monthly means of the above fields for the five years were computed and 
subjected to spherical harmonic analysis. 

It was found from the comparison of the spectral components of 
both sets, M and C, that the climate model generated reasonable 500 mb 
geopotential heights. 

The model temperature field at 850 mb exhibited a generally correct 
structure. However, the meridional temperature gradient was overesti- 
mated and overheating of the continents was observed in summer. 

The structure of the sea level pressure field was not correctly 
reproduced by the model. The four basic pressure systems of the 
Northern Hemisphere were found to be slightly misplaced, with incorrect 
intensitites, and the mass distribution between the hemispheres was 
poorly simulated. 
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Introduction 


During the past 20 years a number of global atmospheric general cir- 
culation models have been developed with the objective of simulating the 
large scale dynamics and physics of the atmosphere. In some cases the 
models were initialized with observed data. and. allowed to run, generat- 
ing a forecast meteorological history. Runs for long time periods 
(years) may be used to compute model- generated climatologies. 

In these experiments tests have been carried out to investigate how 
well- the model simulates the real atmosphere. These tests may consist 
of diagnostics, such as energetics , transports, momenta, and vertical- 
meridional cross sections of various zonal mean atmospheric quantities. 
Forecast statistics and quantitative scalar measures of pattern agree- 
ment, such as correlation coefficients, root mean square errors, and 
gradient "skill- scores," provide additional information for verifica- 
tion, but they do not give detailed information on the quality of a 
model f s simulations of nature. Another verification of model results 
can be done by comparing various horizontal model-generated synoptic 
patterns with the corresponding observed ones through the subjective 
inspection of maps. 

Further information can be provided by comparing the spectral 
components of the above mentioned synoptic fields. These spectral 
components can be represented by the expansion coefficients of a series 
of certain orthogonal functions. In this study these functions will be 
the "surface spherical harmonics." 

Spherical harmonic analysis as a quantitative and objective method 
of analyzing horizontal fields has been applied in the past to many 
geophysical problems (e.g. Spar, 1950; Chapman and Bartels, 1960; Leith, 
1974; Blackmon, 1976; North, 1979). However, spherical harmonic analy- 
sis has not been widely adopted as a method for synoptic pattern veri- 
fication. 


- 2 - 



Spherical harmonics will be used to analyze various horizontal 
fields generated by a general circulation. model (GCM) , the GISS "climate 
model" (Hansen et al., 1979), and to compare the model results with 
nature by analyzing the corresponding observed climatological horizontal 
fields in spherical harmonics. The comparison will be carried- out in 
terms of the relative magnitudes of the spectral components for the 
model and for nature. 

The model-generated horizontal fields are computed by initializing 
the model with global NMC^ data, derived from operational analyses, for 
the first day of a month (at 00. GMT), and, for monthly simulation ex- 
periments, running the model until the end of a particular month, then 
computing the monthly means of various predicted horizontal fields. For 
climatology experiments the model is run continuously for a number of 
years, and the monthly mean fields of the predicted quantities are com- 
puted for all those years. 

By expanding the above horizontal fields in finite series of spheri- 
cal harmonics, the result will be a set of expansion coefficients or a 
set of magnitudes and phases of the spectral components- for selected 
wave numbers, which may be shown in the form of vector diagrams (har- 
monic dials) . 
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Spherical Harmonic Functions 


The spherical harmonics are a complete orthonormal set of functions 
which are solutions of the. Laplace equation in spherical coordinates, 
and satisfy the relationship 



(<j>,A) 
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where 
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ff, is the solid angle which looks at the spherical surface , S , where the 
spherical harmonic functions (SHF's) are defined (dfl = cos<j>d<j>*dA where 
<p is the latitude |“ F > F~j and A. is the longitude J\j, 2 ttJ ) , m is the 
zonal (longitudinal) wave number of the n,m th harmonic, and n - ra is 
the number of nodal parallels, or twice the meridional wave number. 

Both n and m are integers with c£.m^n, while superscripts £ , t in (1) 
indicate the evenness or oddness of the spherical harmonic function 
under the substitution A"* - A . 

The normalized SHF's (see Blackmon, 1976) can be written as follows: 
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m = 0 


where P^ (sincj) ) are the Legendre polynomials (zonal harmonics), and 



P (sixuj) ) are the Legendre associated functions (tesseral harmonics). 
n,m 

The superscripts (e) and (o) indicate the evenness and oddness of the 
SHF's, the eigenvalues n,m being the same as in (1). The normalization 
of the- SHF's is done in such a way that the mean. square value of the 

Y over the surface of a sphere is l/4tr* (see Appendix A) , 

n ’ m (e) (o) 

As an example, the distribution of Y^ 2 and Y^ 2 over the earth is 

shown in Figure 1. As we see, Y, . is a meridional wave with two nodal 

parallels symmetric about the Equator (n - m = 2) and two complete waves 

in the zonal direction (m = 2) . The distributions of Y with lati- 

tude for n<6 are shown in Figure 2. ' The q are zonally symmetric 

functions of latitude with n nodal parallels, and they are variable 

polynomials of degree n. Generally, if n - m is an even number, the 

Y is always symmetric about the Equator since the functions P for 

n,m n,m 

n - m even (odd) are even (odd) functions of latitude. (For more infor- 
mation, see Hobson (1955) . Also, the practical advantages and limita- 
tions of Legendre polynomials and functions for latitudinal expansions 
as compared with Fourier modified series and Chebyshev polynomials, in 
general atmospheric problems on a sphere, are discussed by Boyd (1978).) 
Any horizontal field, Q, which is distributed on a sphere can be ex- 
panded as a function of latitude <f> and longitude X into finite series 
of SHF's as follows: 
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Y v , Y are defined by (2), N is the truncation degree of 
n,m n,m 


where 

n 3 m ii ? in 

the series, and C , S are the normalized expansion coefficients of 
n,ra n,m 

the series (3) . These coefficients can be calculated from the data, Q, 
by using the orthonormality property (1) of the SHF's. (See Appendix B.) 
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The normalization makes it possible to compare coefficients of different 

order m and degree n in the same series or among different series. 

As we see from (4), the truncation of the series at N = k does not 

affect the calculations of the coefficients C , S for msn = k be- 

n,m n,m 

cause in the above integrals the independent variables are not the n,m 
but 4*> X • 

We can calculate the amplitude A and the phase angle B of 

n,m r n,m 

any harmonic from 


n,m 



2 

+ S 
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The above relations are as indicated in the harmonic dial in Figure 3. 
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Fig. 3. 


For calculating the expansion coefficients C • , S , a closed-- 

n,m n,m 

form numerical integration scheme is required to approximate the in- 
tegrals (4) . Since the data, Q, are given on grid points (discrete 
function of <J> , X ) , this makes necessary the discretization of the Y 

n,m 

at the same grid points, which leads to the evaluation of the Legendre 

and Fourier sine-cosine functions on the above grids. 

Generally tables of normalized associated Legendre polynomials are 

available (e.g. Belousov, 1962) but they do not correspond to the GISS 
o o 

grid spacing (8 x 10 latitude- longitude grid). So, following the in- 
structions of Belousov (1962) and Hobson (1955), the Legendre poly- 
nomials and associated functions were computed by means of Rodrigues* 
formula (6) , 


P n , m (sin<|>) = (cos<j> ) 
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up to n = m = 18, for the appropriate latitudes <f>, and stored on disk 
for use as required. The calculations were done by extended precision 
(28 significant figures) on an IBM 360/95 computer. 

An attempt to recreate Belousov’s tables by using the algorithm (6), 
for the same latitudes’ up to n = m - 18, showed perfect agreement. (Al- 
ternative numerical methods for computing Legendre functions and poly- 
nomials are discussed by Merilees (1973) .) 

Various numerical integration schemes for the calculation of the 
integrals (4) were tested by checking the orthonormality property (1) 
of the SHF's. This was done by substituting Q( «}> , X ) in (4) for one 


of the Y (e) , Y ( °\ 
n,m n,m 

Finally, the method that has been adopted was a combination of the 


Simpson and "three-eighths" rule (Krulov, 1962) for the <p integration 
and Filont's rule (Davis and Rabinowitz, 1967) for the X-integrals. . 
The integration was carried out in two steps: 


TT/2 

a) ^-integration F n m (X) = R n/ia * J Q ($ , X) * p n ^ m (sin<|>) *cos<{>d4> (7) 

— 71/2 
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With the above method, equation’ (1) was satisfied up to the sixth deci- 
mal place. 



Examples and Tests of the Method 


The area-weighted mean value, Q, of any synoptic field, Q, is 


i 

1 ° -h J Q{ *' 

s 

i 

From (2) we have Y - 

v 7 a n 


X) *dfl . 


(4ir)'^; so 


Q = -► 

s 
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dQ = Y »C . Therefore 
O 0,0 0,0 


(9) 


A coefficient C represents the amplitude of a zonally symmetric 

n ’° •*% 
function of latitude only. For example, Y = (4/3. tt ) 2 . sin 43 , 

-3- 2 i5U 

Y - (3.2-rr) 2 . (3-sin4> - 1). Then we may say that the amplitude 

W J v 

of the mean meridional difference between values of Q at the North and 
South Poles is approximately represented by Also, C 0 ^ is a 

rough measurement of the mean meridional Fole-to-Equator Q-field dif- 
ference . 

By using the already described numerical methods, the observed 

climatological fields of temperature at 850 mb (T850) , 500 mb geo- 

2 

potential height (G500), and sea level pressure (SLP) were analyzed 
in spherical harmonics. Each field is represented by a set of 361 
nontrivial expansion coefficients. From the series (3) the fields Q 
were reprod uced for different truncations N, up to N « 18. 

The observed climatological fields were obtained from NCAR, the 
National Center for Atmospheric Research, in Boulder, Colorado. 
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Figures 4, 5, and 6 indicate three climatological fields for the 
month of December which were subjected to harmonic analysis. 

An area-weighted root-mean-square (rms) error test (10) was carried 
out between the initial fields, Q, and the reproduced ones, as fol- 
lows: 





( 10 ) 


In Figure 7 a, b, c are plotted the rms errors of the reproduced December 

o 

climatology for T850 in degrees G , G500 in meters, and SLP in millibars, 

3 

as a function of the truncation N. As we see in Figure 7 a, b, there 
is a big reduction of the rms for 1 < N $ 2 due to the symmetrical zonal 
structure of T850 and particularly G500. By truncating the series at 
N - 2, the hemispheric and the zonal structure (Y_ , Y„ ) of those 

* ) v ^ j U 

fields has been captured* Therefore, the Y 9 harmonic accounts for a 
large percentage of the total variance of T850 and G500. The small 
scale features of the three fields can be captured with truncations 
N-^8. We observe that for these truncations (Figure 7) the rms error 
does not diminish at the same rate as for This leads to the con- 

clusion that only a small portion of the variance in these fields can 
be explained by using harmonics up to n = m - 18. 

The existence of a non- zero rms error, for N = 18 is characteristic 
of the integration "noise." The noise can be reduced by reproducing 
the fields, Q, from only the leading harmonics. This was done by rank- 

ing the co efficients in descending order and reproducing the field using 
3 

The number of harmonics which is used each time for the reproduction 

2 
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only the first 10, 20, 40, 80 dominant harmonics. As before, an rms 
error test was carried out and the results are illustrated .in Figure 8 
a, b, c. 

For the above reasons it will be instructive in the course of this 
study to represent the fields by a small set of coefficients, corres- 
ponding to their leading harmonics. 

Figures 9 a', b, c, d - 11 a, b, c, d show how the fields are 
simulated when reproduced from 10, 20, or 40 leading harmonics, and 
all the harmonics. Note that by adding more harmonics, small wave 
length disturbances become apparent in the synoptic maps. 
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The Model Experiment 


A global general circulation model has been developed for cl ima te 

circulation studies by a group under J. Hansen at the Goddard Institute 

for Space Studies (GISS) in New York. This model can be operated with 

o o 

a horizontal grid of 8 of latitude, <p , by 10 of longitude, X , and 
seven sigma ( cr) layers in the vertical. With this resolution one 
simulated month takes about two hours on an IBM 360/95 computer. (See 
Hansen et al., 1979 for further details of the climate model and Somer- 
ville et al., 1974 for a description of- the first GISS model from which 
it was developed.) 

In order to generate a climatology history, the GISS climate model 
was initialized with observed global data derived from operational NMC 
analyses and interpolated to the GISS grid. The initial conditions 
were for the first day of December 1976 at 00. GMT. Then the model was 
allowed to run for five simulated years until the end of January 1982 
and the monthly outputs were averaged to obtain a five-year model 
climatology. 

The model- generated (M) global fields of sea level pressure (SLP), 
geopotential height at 500 mb (G500) and 850 mb temperature were se- 
lected for verification. The monthly means of these fields for all five 
years were analyzed in series of spherical harmonics. 

Global observed climatological data (C) of the same horizontal 
fields provided by NCAR and interpolated to the GISS grid were also 
subjected to spherical harmonic analysis. Each field was represented 
this way by 361 expansion coefficients for the "global" harmonic analy- 
sis. 

In addition, by assuming that the data, M and C, are even functions 
of latitude about the Equator (reflecting the Northern Hemisphere data 
into the Southern Hemisphere), we can eliminate possibly incorrect data 
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due to the lack of observations in the Southern Hemisphere, and apply 
the spherical harmonic analysis by' calculating the integrals (4) and 
the relations (5). This time, since we have assumed that the data are 
even about the Equator, only coefficients with even n - m will be non- 
zero (see Appendix B) , and each field will be represented by 190 expan- 
sion coefficients. This will be called a Northern Hemisphere analysis. 
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Results 


In order to be able to compare the model-generated climatology (M) 
with the actual climatology (C) , it was assumed that the actual clima- 
tology remains in steady state for the next five years. 

Interpretation of the results will be done: 
a) by tabulating and ranking the leading harmonics of each field in 
descending order of the magnitudes of the amplitude coefficients. 


A ; 

n,m 

b) for a particular field, by comparing the annual cycle of each 
component (amplitude and phase of the harmonic) for M with the 
cycle of the same component for C. 

The annual cycles of only the leading components for the three, 
synoptic fields will be examined. Also, seasonal climatologies for M 
and, C were obtained by averaging the amplitudes of the harmonics over a 
three-month period (see Appendix B) . The seasons were defined as fol- 
lows: 

Spring : March, April, May; Summer : June, July, August; 

Fall : September, October, November; Winter : December, January, February 

For the above seasons the eight leading harmonics with their magnitudes , 

A , and their phases, p , will be tabulated. For the case where 
n,m n,m 

m = o, the phase angle p is defined as follows: 

n,o 


180° 

if 

c < o • 

n, o 

0° 

if 

C >0 

Xl/O 


(ID 


The synoptic fields may then be represented in terms of the amplitudes 
and phases of the dominant components of each field. 

Table 1 contains the eight leading harmonics, not including A n n . 
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■ Table 1 


Global mean valxies ahd the 8 largest 'hatmbnics A for each field 

71 / m 

(T850, G500, SLP ) for winter actual climatology (C) and winter model 

generated climatology (M) 

(global sph. harm, analysis) Code: .n,m 

A • / 0 
n,nr *n,m 


785 0 

M EA N 

1 

2 

3 

4 

5 

6 

7 

* 

8 

\ * 


- 2*0 

1 * 0 

3*2 

. 2.1 

3 * 1 

4*2 

5,2 

5* 4 

C 

279*7 

40# 1/180 

1 1# 4/1 80 

3*6/ 73 

3.4/341 

3*1/306 

3*0/ 53 

ro 

*• 

CO 

* 

CVI 

2*6/333 



.. 2*0 

1 f 0 

6.0 

3*1 

1 3 *2 

2*1 

8*0 

6* 2 

’ M 

warn 

4 3 *8/ 1'8 0 

• 12*1/130 

5. 2/ 0 

4. 1/302 

3.9/ 77 

3*5/340 

3. 4/180 

2.8/ 43 

* 

csoo 

MEAN 

1 

• 

2 

3 

, 4 

5 

' a ‘ , 

7 

8 

, 


.2*0 

6. 0 

' 4*0 

1 .0 

3.0 

7,0 

2,1 . 

4*1 . 

c 

5636*0 

9 16*4/180 

112.1/ 0 

11 1 . 0/180 

83.9/130 

76.4/ 0 

59. 5/180 

51.9/ 10 

49*5/ 35 


• . 

2*0 

1 » 0 ' 

6*0 

4.0 

4,1 

'.2,1 

3,1 

5 * 0 - 

M 

5605*2 

843.0/180 

191. 2/180 

94.3/ 0 

85.3/1 80 

53.0/ 10 

47* 0/358 

46. 5/336 

42*2/ ,0 

* 

SLP 

MEAN 

1, 

2 

3 

4 

5 

6 

7 

a 

‘ - 


1 *0 

2 9 0 

4.0 

6,0 

’ 3*0 

7,0 

4,2 

2. 1 

* 

- c 

1011 .2 

18.5/ 0 

1 C* 7/1 80 

1 0. 3/1 80 

8.5/ 0 

7.2/ 0 

4. 8/130 

4.7/192 

4.4/ 88. 



A *0 

1*0 

2*0 

3, 1 

8 * 0 

3, 2 

2,1 

4* l’ 

' M 

1010*7 

9 . 4/1 80 

9*3/ 0 

.6*6/ 0 

4 . 9/ 1 09 

4*6/ 0 

4..6/241 

4.5/ 97 

3.9/ 65 












































for the winter C and M fields. Included in Table 1 are the harmonic 

indices, n,m, and amplitudes, A ', and the phase angles for the eight 

n,m 

leading harmonics. Also tabulated are the global mean values of each 

field, computed from A using equation (9) . Reasonable global means 

o 5 o 

•were generated by the model for T850, SLP and G500. The same dominant 
•harmonics, 2,0, for G500 and T850 appear in M and C. However, the model 
has overestimated the polar equatorial difference, 2,0, in T850 and 
underestimated it in G500 and SLP, with a phase difference of 180° in 
.the last field. The difference between the poles, 1,0, is grossly 
underestimated by the model in SLP (18.5 vs 9.3) and overestimated in 
G500. The 1,0 in the case of SLP indicates an excess of mass in the 
Northern Hemisphere during winter-. 

The differences between the eight leading harmonics for C and M 
indicate the character of the predicted climatology error which is' 

dominated by the 6,0 harmonic in T850 (not found in the C list) and in 

SLP (not found in the M list). Also in SLP the error is in the 3,0 and 
7,0 terms (contributions to the asymmetry of the two hemispheres) which 
do not appear among the eight leading M harmonics. 

The data listed in Table 2 are for the spring. Again the global 

means for T850, SLP and G500 were correctly estimated by the model. 

The polar equatorial difference is still dominant and 4,0 is the second 
dominant harmonic for T850 and G500. The 2,0 is again overestimated by 
the model in T850 but is well estimated in G500. Also the 4,0 is well 
estimated for T850 and G500. In SLP the 2,0 is overestimated by the 
model with a phase difference of 180 . Also overestimated are the 7,0 
and 8,0 terms. Additionally, the 1,0 and 3,0, which are the dominant 
harmonics in C, do not appear among the .eight leading harmonics in M. 
Generally the M error in SLP is dominated by the 5,0 and by 1,0, 3,0, 

4 jO, which are not found in the M list. 

i 

Table 3 contains information as in Table 2 but for the summer 
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Table 2 > 


Global mean values* and ‘the 8 largest harmonics ’A m for each field 
(T850 / G500 / SLP) for spring actual climatology (C) and spring model 

generated climatology (M) 

{global sph. harm, analysis) Code. n,m 


” n,m 


A / B 
n,rtr M n,m 


7850 

' MEAN 

C 

280 .7 

M 

280.0 

G50 0 

MEAN 

C 

5651 .2 

M 

5618.6 

I SL p 

MEAN 



ion .0 


ion «o 



2,0 

40*5/ 180 
2 » 0 

46. 1/180 


4,0 

4# 7/1 80 
l 4.0 
4. 8/1 80 


2,r 

3*5/ 42 
5,0 

4, 3/ 0 


5.0 

3.4/ 0 

2.1 

4.1/ 34 


,1 ,1 

2.8/ 2*7 
8,0 

3.7/180 


5,4 

2. 3/340 
7,0 

3. 6/ 0 


8,0 

2.0/180 

1,1 

3*5/ 39 


3, 2 

1.9/127 

6,2 

3.3/ '20 



889.0/180 122.2/180 123.3/ 0 119.1/ 0 77.0/ 0 4 4. 9/ 31 I 37-.6Z 0 137.5/ 38 
2.0 4,0 6,0 6,2 7,2 . 5,1 I 5,2 I 5,0 

.33U8/180 1 3 7o 6/1 80 66.2/ 0 37.2/ 50 36.2/ 65l 33.3/323 1 33.2/ 61 133.0/ 6 



15.8/ 0 1 12.5/ 01 10.2/130 


2,0 

12.4/ 0 


8,0 

7.3/ 0 


7,0 

7.5/180 


4.0 

9.3/180 

5.0 

6. 0/130 


6*0 


77.0/ 

0 

7,2 


36.3/ 

65 

5 

6 ,0 


8.6/ 

0 

8 * I 


. 5*1/255 



7,0 

4* 6/ 1 8 0 
6 , 0 

5, 1/ 0 


‘ 8,0 
3.9/ 0 

'5,1 

4.6/ 79 


8 , 1 

3.5/244 

1 , 1 

4. 1/230 












































Table 3*. Global mean values and the 8 largest harmonics A for each field 

n,m # 

(T85Q, G500, SLP) for summer actual climatology (C) and summer_ model 

generated climatology (M) 

(global sph. harm, analysis) Code. n,m 



£82*4 


1010.6 ' 


• 2 , 0 I 1,0 I 4*0 
32. 4/1801 21.2/ oj 6,3/180 


3S. 0/1801 25.6/ 0| 10.0/180 



2 

• 1,0 


4 7€. 5/ 

'o 

1 , Q- 


462.4/ 

0 

, 



2 » 1 

4.3/ 73 
_ 2 » 1 
6.3/ 68 




14. 5/ 180 I 14.2/ 


3 

l »Q 


1 0. 0/ 

0 

, 6,0 


8.6/ 

0 



6*0 

8. 1/ O 
, 2.0 
3.3/ 0 


5 

6 

to 

5,0 

2.6/218 

2. 6/ O 

7,0 

4,2 

4,7/ 0 

,4, 7/166 

5 

6 

6*0 

9,0 

74.4/ 0 

37, 1/130 

3,0 • 

7,0 

46.3/ 0 

33.0/ 0 

5 ' 

6 

4*0 

i.i 

8,1/180 

6.3/293 

5,1 

5,0 , 

6.9/ 72 

6* 8/1 30 



7,0 

2.3/ 0 

* 

3,2 

4. 7/168 


2,2 • 
2. 1/143 
1, 1 

3*4/ S5 




2,1 

6. 1X262 
3,2 

6 , 8 / 1 


4,1 

4.4/ 61 

1,0 

6.7/180 
















































period. The global means were very well estimated by the model for the 
three fields. Again the dominant harmonic for T850 and G500 is the 2,0, 
due to the strong meridional gradient of these fields. In G500 the am- 
plitude o'f this gradient is underestimated by the model. The other 
-dominant components for the G500 and T850 are the 1,0 and 4,0. All of 
them except the 4,0 in G500, and especially in T850, are overestimated. 

From winter' to summer we find, comparing Tables 1 and 3, an in- 
crease in the relative magnitude of the 4,0 harmonic, notably in T850. 
This appears to reflect the heating of the mid-latitude continents in 
the Northern Hemisphere summer, as a result of which high temperature 
and geopotential values are observed north of the geographical Equator. 

This meridional pattern is represented by the zonal harmonic 4,0 with 

o 

the phase angle- of 180 , indicating relatively high temperatures and 
geopotential in middle latitudes. The 3,0 term which contributes to 
the asymmetry of G500 in both hemispheres is badly estimated by the 
model. The M error in this case is dominated by the 4,2 in T850 and 
by 6,0 and 5,0 terms in G500. In SLP the 2,0 harmonic is significantly 
•underestimated with opposite phase as is the 1,0 term. The M error is • 
dominated by the asymmetric term 3,0 and by 4,0. The other harmonics, 
6,0 and 2,1, are found in both the C and M sets with about the same 
amplitudes and' phases. 

Fall results are indicated in Table 4, which shows that nearly 
perfect global mean values were estimated by the model for all fields. 
The 2,0 term is dominant in T850 (overestimated) and G500 (underesti- 
mated). The 1,0 and 4,0 are the next dominant components in T850 and 
G500. Beyond the third dominant harmonic 4,0 there is little corres- 
pondence in T850 between the C and M harmonics.' In SLP the 1,0, 2,0 
and 3,0 terms are grossly underestimated with opposite phase in 2,0. 

The 4,0 harmonic is not found in the M list while the 7,0 and 5,0 are 
not found in the C list. 
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Table 


Global mean values and the 8 largest harmonics A n m for each field 

(T850, G500, SLP) for fall actual climatology (C) and fall model 

generated climatology (M) 

(global sph. harm, analysis) Code. n f m 

^n,ii/ ^n,m 



5631 • 3 


804.4/1301 205.3/ 01133. 7/160 153.2/ 


33.0/ 0 3 1.9/ 83 30.2/305 1 29.7/ 29 


SLP 

MEAN 

1 

• 


1,0 

c 

1011.0 

15.4/ 0 



2,0 

M 

1010.7 

9.4/ 0 


'Z 3 4 5 

2 » 0 3,0 4,0 6,0 

7/130 12. 5/ 0 12. 4/130 7.6/ 

1,0 7,0 8,0 3.0 

5/ 0 .5. 8/180 5.7/ 0 5.4/ 





















In the spherical harmonic analysis for the Northern Hemisphere 
alone, the asymmetries between the two hemispheres are. eliminated and 
the effects of oceans and continents in T850 and SLP are more pronounced. 
The ocean- continent distribution in the Northern Hemisphere corresponds 
to zonal wavenumber m - 2. For example, -during the Northern Hemisphere 
winter the oceans are warmer than the continents (cont mentality) and 
zonally the temperature distribution consists of two alternating high 
and low values due to the two large continental areas (North America, 
Eurasia) and two oceans (Atlantic, Pacific). The above pattern corres- 
ponds to the harmonic 4,2 as indicated in Figure 1. The same pattern, 
with a J small phase shift, fits the SLP field due to the longitudinal 
distribution of the four basic pressure cells in the Northern Hemisphere 
(Aleutian low, Canadian high, Icelandic low, and Siberian high). 

Table 5 contains the leading harmonics of the three Northern Hemi- 
sphere' fields for winter. The model has generated a reasonable hemi- 
sphere mean value, for T850, with a slight cold bias. The area-weighted 
means of G500 and- SLP- are underestimated. Differences between the North 
Pole and Equator, 2,0, are overestimated by the model, especially in 
SLP, for which 2,0 is not even found among the eight dominant harmonics 
.in the observed climatology (C) . Although the model generates a 2,0 
harmonic as the leading component of SLP, in nature this term does not 
even appear among the eight leading harmonics. The 4,2 is the dominant 
SLP harmonic in C and the second dominant terra in T850. The model 
generates a reasonable 4,2 harmonic in SLP and-T850, although it is not 
the leading terra in M. This means that the locations and the intensi- 
ties of the above pressure cells were not correctly simulated by the 
model. Also, the 4,2 in T850 was underestimated. The M error in the 
T850 is dominated by the 6,0. The harmonics 6,0, 3,1 appear among 
the dominant terms in the M G500 but in different order. The 7,1 and 
6,0 represent model climatology errors- in SLP. 
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Global mean values and the 8 largest harmonics A for each field 

xi f m 

(T850, G500, SLP) for winter actual climatology (C) and winter model 

generated climatology (M) 

(Northern Hemisphere analysis) Code: ' n,m 


' n,m 


5 

6 

6,2 

2,2 • 

O 

< 

O 

• 

3.8/ 68 

4,2 

6, 2 

5.6/ 55 • 

5.0/ 31 

‘5 

6 

4,2 

6,2 

74*7/ 96 

6 4.3/ 6 9 

5,1 

5,3 

77.1/ 14 

64. 7/283 

5 

6 

6,0 

5,1 

5.7/ 0 

4. 9/177 

2 1 a 

6,2 

5.4/222 

5.4/173 




5*3 

3.3/276 
1 ,1 

4.7/ 2 


5. 1 

2.7/ 58 
7,1 

4.4/ 70 



1,1 - 


5612.-4 I 919.6/180 1 85.8/ 0 182.2/358 78.9/315 74*7/ 96 |64.3/ 69 I 58.5/180 |53*8/339 


948.5/180 llQO.O/351 93.4/ 0 |78.5/ 75 177.1/ 14 64. 7/283 I 55.7/ 79 154.1/323 





1012.7 


4 ♦ 2 

8.7/193 

2,0 

13.4/ 0 


4, 0 

7. 4/180 

3,1 

9.8/1 00 


, 3.1 

7. 2/106 
4,0 

8. 3Vl80 


7.1 

6. 3/226 

4.2 

7. 2/207 


• 6,2 
4. 7/186 

1 ,1 

3. 7/122 


2 , 2 

4.7/208 

7,1 

3.6/267 













































• Table 6 contains the data for the spring. The same dominant har- 
monic 2,0 appears for all fields, and it is grossly overestimated by 
the model, especially in SLP. All the area-weighted means are under- 
estimated as well as the majority of dominant harmonics in T850. The M 
T850 error is dominated as before by the 6,0 and 3,1 terms. In the 
G500 the 4,0 and 6,0 are overestimated. In SLP the 7,1 and 6,0 are 
underestimated as in winter. The 6,0 is not found among the first 
eight harmonics in M. This is considered to be an M structure error 
in SLP because the 6,0 corresponds to the C meridional pressure dis- 
tribution shown in Figure 12. 



Figure 12 : Meridional pressure pattern at sea level 

(From J. Spar, Earth, Sea and Air , 1962) 
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Table g. 


Global mean values and the 8 largest harmonics A for each field 

n,m 

(T850, G500, SLP) for spring actual climatology (C) and spring model 

generated climatology (M) * 

(Morthern Hemisphere analysis) Code; n,m 

A n,n/ ^n,m 


785 0 

MEAN 

C 

M 

2 81 .1 
2 80*2 

QS0 0 

MEAN 


5671 .1 



2 ,0 

39.4/180 

2,0 

47.4/180 


. 1,1 
6.1/ 36 
6 , 0 

£. 1 / 0 


2,2 

.3. 2/1 25 

1 , f 

7. 1/ 41 



6,2 

2.8/ 28 
6,2 

6 . 6 / 13 



752.0/180 I 70.7/ 0 152. 4/ 


827.4/180 1101.9/1 80 89. 0/ 



10 10.7 


2,0 

6.1/ O 

2, 0 

14.7/ O 


7, 1 

£.4/241 

1*1 

7. 6/240 


6,2 

4. 6/177 
6,2 

6*4/1 63 


6,0 

4.6/ -0 

5,1 

5. 2/ 34 


5 

6 

4 t 4 

4,2 

2.2/352 

2. 0/1 01 

2.2* 

3,1 

4*3/120 

4. 2/334 

5 

6 

4 ,0 

8.0 

4'7. 6/180 

46. 5/180 

3, 1 

7,3 

47.0/351 

4 3. 9/336 

'5 

6 

i , i 

4.0 

4.4/257 

3. 2/180 

4,0 

2.2 

4.4/130 

4. 1/3 06 



8,0 

1 . 9/180 ’ 
5,1 

3. 8/285 


9. 1 

1 .9/2 17 
4, 2 

3.7/ 93 




8,4 

2.4/ 28 
3,1 

3. 8/143 


2.4/272 
7, 1 

3.4/247 














































For the summer (Table 7) , the mean T850 is warm-biased and there 
is a deficit of mass in the Northern Hemisphere, indicated by the under- 
estimation of the mean SLP, while the mean G500 is well simulated. 

The dominant harmonics for T850 and G500 are the 2,0, which were 
correctly simulated by the model, and the 4,0, due to the insolation, 
which is grossly overestimated in both M fields. _ The same tendency of 

4.0 in the summer was observed in the global analysis. 

The 4,2 is overestimated by the model in T850 and SLP with a phase 
shift in SLP. This indicates that during summer there is an opposite 
continentality effect which also results in an opposite zonal pressure 
distribution than the one in the winter. The above is clearly indicated 
by the increase of 4,2 in SLP, with 180° phase shift from winter to sum- 
mer. A monsoonal pressure distribution term, 1,1, appears in both' M and 
C as the leading harmonic for SLP. The 3,1 and the 2,2 which represent, 
together with 4,2, the continentality effects are overestimated by the 
model in both T850 and SLP. The M error in SLP is dominated by 2,0 
(not in the M list) . 

For fall (Table 8), we have reasonable means for T850 and G500. 

The same components, 2,0 and 4,0, as in the summer are dominant in the 
above fields. In T850 there is an overestimation of the 2,0, 4,0 and 
6,0, while the 1,1 and 3,1 are correctly estimated. The M error is 
dominated by 6,2. 

In G500 the 2,0, 4,0, and 6,0 terms are very well simulated by the 
model, while the 6,2 and 4,2 harmonics are overestimated in SLP. The 

2.0 dominates in the M SLP and represents the largest "model climatology" 
error, because the above harmonic is not found in the C list. Also, the 

7.1 is not among the leading harmonics in the M list, while 4,0, which 
dominates in C, is underestimated as is 3,1. 
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Global mean values and the 8 largest harmonics A for each field 

n ,m 

(T850, G500, SLP ) for summer actual climatology (C) and Rummer model 

generated climatology (M) 

(Northern Hemisphere analysis) . Code: n,m 


■ 785 0 

MEAN 

C ‘ 

2 87.4 


2 89.1 ‘ 

G50, 0 4 j 

MEAN 

1 * 




* 

5 

6 

3.3 

•2,2 

3.2/ 87 

3. 1/152 

2,2 

3. 1 

' 6.4/156 

5. 2/ 7 3 

5 

6 

* 

10 » o 

5,1 

32.6/180 

3 0. 9/132 

3,1 

7,5 

1 8.7/114 

1 8. 1/241 

> 


5 

6 

' 6,2 

7.1 

4.5/149 

3. 7/236 

3,1 

5,1 

7.0/247 

6.8/ 63 


5782.2 


5769.1 


• M54N 


1011 .6 


10 03.6 




20 .8/ 180 1 e.5/1 80 


2 2.4/180 16.2/1 80 


1.1 

5.4/ 54 
4,2 

8. 6/173 


5,1 

3. 2/208 

1 , 1 

• 8.2/ 63 


6 , 0 - 
3.0/ 0 

•5,1 • 
4. 5/222 



.4*2 
3.0/173 
8 . 2 

4.2/359 





•/ 1,1 2 , 2 . 
11.7/275 €.3/356 


1 2.8/255 11.7/ 6 


5, 1 

5.9/ 80 
2,2 

7. 9/350 


4,2 

4.8/ 38 
2 , 0 . 

7.9/ 0 


3,1 

2. 9/232 
• 8.2 
4.3/183 


4,0 

2.6/180 

6,4 

4o 1/196 
















































Global mean values and the 8 largest harmonics A^ M for each field 

, . n,m 

<T850 , G500, SLP) for ,fall actual climatology (C) and fall model 

•generated climatology (M) 

. , - . , , - r . Code: . n,m 

(Northern Hemxsphere analysis) 

A / 6 

• n,nr M n,m 


MEAN 


2S2 *5 


2 82*5 


MEAN 


M EAN ' 


2,0 

1 

34*8/180 
2,0 .• 
40.3/ 180 


4,0 

4* 2/1 SO 
' 4, 0 

e'. 3/.1 eo 


i,i' 

4* 1/ 12 

6,0 

a. 5/ 0 


3* 0/ 0 

6.2 

4.5/ 12 


3 , 1 

3.0/339 

1 . 1 

4.2/ 48 


4,2 

2.5/ 85 

3,1 

2. 9/345 


6,4 

2. i/302 
.5,1 * . 
2* 6/299 


7,1 

2* 1/ 54 

7,1 

2 . 6/ ' 48 


1 

i 

2 

3 

4 

5 

6 

7 . 

8 

2,0 

699. 2/180 

4,0 

153.1/1 80 

6*0 

88» 9/ 0 

‘ 3* 1 ‘ 

63. 4/ 6 

6,2 

4 6.0/ 95 

4*2 

34* 6/9 7 

r,i 

33. 0/352 

7,3 

28*8/ 18 

2,0 

690.3/ 180 

• ■■■■■■■■ ■ 

4,0 

176.2/180 

/6 ♦ 0 

79* 7/ 0 

6 , 2 

.56.5/ 60 

7,3 

43.6/348 

4*2 

4 1*6/ 54 

5,1 

35*3/ 14 

8,2 

32*6/ 75 




4 .0 

7*7/180 

2.0 

17.5/ 0 


7, 1 • 
£.0/233 
‘ 5, 1 

4.4/ S3 


' 3,1 
4. 1/1 01 
• 6,2 
4. 3/165 


6,2 

3. 7/ .173 
•3,1 
3.6/145 


6,0 

3.1/ 6 

4 » 0 

3.3/180 


* 8 , 0 
2 . 9 / 0 

1,1 

3. 2/241 


4,2 

2* 3/186 
"6,4 
2. 8/177 


5,1 . 

2.6/125 

6,0 

2.7/ 180 

















Figures 13 and 14 illustrate the annual cycles of the components 
A q ^ of T850 , G500, and SLP for the global and hemispheric analyses 
• respectively. The dashed curves represent the model-generated spectral 
components and the continuous curves correspond to the actual climato- 
logical data. The global means are related to A through relations 

o >o 

(9) and (5). Therefore, the annual cycle of A q in Figures 13 and 14 

i$ proportional to the annual cycle of the global means for the above 

fields. The units for the amplitudes A are degrees Celsius for 

n,m 

T850, meters for G500, and millibars for SLP. To find the global or 

hemispherical mean values, the plotted A values must be multiplied 

V 4 0,0 

by (4-tf )” 2 or 0.28. 

In Figure 13 there is an annual cycle of the global mean C T850, 
with maximum in July and minimum in January, due to the large percentage 
of land in the Northern Hemisphere, which dominates the global annual 
temperature cycle. The climate model generates a similar annual cycle 
for the mean T850, but with a larger amplitude (overestimate in summer, 
underestimate in winter) . This can be interpreted as a tendency of the 
model to overheat the continents in summer and overcool them in winter. 
This is also indicated by the annual cycle of 0,0 for T850 in the 
Northern Hemisphere (Figure 14) , which, as expected, shows a larger 
amplitude than the global mean for both C and M. A similar annual cycle 
is observed for the mean G500. The model underestimates the winter mean 
values of G500 while it successfully captures the summer means. 

There is virtually no annual cycle of the observed global mean SLP 
(Figure 13) . This is correctly represented by the model with only a 

small diffe rence in the winter."* However, the distribution of mass 

4 

For A^ the units of SLP, G500 are mb - 1000, m - 5000, respectively. 

5 

The model is artificially corrected during the climate run to maintain 
constant global mass. 
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Figure 13. 
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between the Northern and Southern Hemisphere is not correctly simulated 

by the model. This is indicated by the existence of a biased annual 

cycle in the mean SLP for the Northern Hemisphere (Figure 14) . In the 

observed climatological SLP fields, the distribution of mass between 

the hemispheres is such that always the mean pressure is greater in the 

Northern Hemisphere. At least in winter this is associated with the 

large percentage' of land in northern latitudes, which results in the 

generation of continental high pressure systems. The distribution of 

mass between the two hemispheres is indicated in Figure 15, where the 

1,0 spectral component is plotted. The climatological C 1,0 harmonic 

shows that the Northern Hemisphere has maximum of mass in January and a 

minimum in July, with an excess of mass over the Northern Hemisphere all 

year, as shown by the constant phase angle P =0. The model simu- 

1, u 

lates the same minima and maxima in the 1,0 harmonic. However, from 

April through August the mean pressure in the model climatology is 

greater in the Southern Hemisphere, as indicated by the change of phase 
o o 

angle from 0 to 180 . This appears to be associated with excessively 
high pressures in the Antarctic in the Southern Hemisphere winter. 

The annual cycle of 1,0 for T850 (Figure 15), showing the seasonal 
temperature difference between the hemispheres reflected in the phase 
shift, is rather well simulated by the model but with a small phase 
delay in April and a temperature excess in the Northern Hemisphere in 
July and August, indicated by the large amplitude of the 1,0 harmonic 
in M. The same is observed for G500, with high geopotential values 
over the Northern Hemisphere in summer, which is correctly represented 
by the model. 

Figures 16 and 17 represent the mean polar-equatorial differences, 
2,0, for the globe and the Northern Hemisphere, respectively. For T850 
and G500 the 2,0 harmonic is related to the distribution of insolation. 
In the 2,0 C global cycle (Figure 16) of T850, maximum and minimum 
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Figure 15, 
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Figure 17. 
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amplitudes are found in March and July, respectively. The M cycle of 
2,0 for T850 follows this pattern but with a high bias, including an 
overestimate of the maximum in March and a second maximum in November, 
indicating larger temperature differences between Poles and Equator. 
Similar cycles are observed in G500 but with a low bias in the global M 
cycle . 

The C -2,0 annual cycle for the Northern Hemisphere (Figure 17) for 
T850 and G500 has one minimum (in July) and one maximum (in January) . 

The model simulates these extrema with a lag of one month for both 
fields. As before, there is an overestimation of the maxima in winter, 
which results in too strong a meridional temperature gradient. In G500 
the cycle of 2,0 for the Northern Hemisphere is very well estimated by 
the model. On the other hand, in the model simulation of the 2,0 har- 
monic of global SLP (Figure 17), the SLP cycle of 2,0 is correct in 
phase but positively biased in amplitude, notably in October. Appar- 
ently the Southern Hemisphere is mainly responsible for the global 
phase error. 

The harmonics 3,0 and 5,0 (Figures 18 and 19) represent, like 1,0, 
the asymmetry of the hemispheres. The annual cycles of these harmonics 
for T850 are poorly represented by the model, with large phase and am- 
plitude errors. For G500 the model annual cycle of 3,0 parallels that 
of the observed climatology but underestimates the amplitudes, while 
the 5,0 is reasonably well estimated by the model. In the C SLP there 
is a significant annual cycle of the 3,0 harmonic, which is not cor- 
rectly represented by the model, indicating an inadequate exchange of 
mass between the hemispheres. 

The annual cycle of the 4,0 harmonic is plotted in Figures 20 and 
21. As mentioned before, the 4,0 term for T850 is related, like the 2,0, 
to the distribution of the insolation. For the global C cycle we have 
one minimum in January and one maximum in August. The latter lags 
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Figure 21, 


TIME 'SPECTRA FOR THE COMPONENT 4,0 

TB50 

AMPLITUDES A . : PHASES B 


JFMAMUJASOND 


JFMAMJJASOND 


AMPLITUDES A . 


PHASES B . 


U F M A M J J 

A S O N D \ 


'JFMAMJJASOND 



SLP 


AMPLITUDES 

A . 


PHASES B . 


JFMAMU.JASOND 


JFMAMJJASOND 


HKTTEFM IBlISRifFE; ANALYSIS 


tf * * OBSERVED CLfMATCLOGr 

0 0 * PREDICTED CLIrtATCLCGT 












behind the June solstice due to the thermal inertia of the earth's at- 
mosphere. The model displays the same annual cycle, but the larger am- 
plitude due to overheating of the continents in northern latitudes. 

Also, the maximum in the M T850 cycle occurs in July, -which indicates 
too little thermal inertia in the model atmosphere. A similar cycle 
for 4,0 in T850 is observed for the Northern Hemisphere (Figure 21), 
with two maxima and two minima associated with phase shifts. This 
cycle is correlated with the annual insolation. The winter cycle of 
4,0 for T850 is perfectly simulated by the model, while the summer 
cycle is overestimated. This result further illustrates the fact that 
the model artificially overheats the Northern Hemisphere in summer. 

This is also indicated by the large amplitudes of the model (4,0) of 
G500 in summer (Figure 21) . For SLP, the amplitude of the global 4,0 
harmonic (Figure 20) is underestimated, with incorrect phase from May 
through September, while in the Northern Hemisphere (Figure 21), the 
model more closely parallels the observed amplitude, although the phases 
are still shifted in summer. 

The 6,0 harmonic (Figures 22 and ,23) is important as it represents 
the SLP meridional structure. However, in both the global and hemi- 
spheric analyses, the M annual cycle of 6,0 for SLP is very different 
from the corresponding C cycle. This may be due to the errors found 
over cold areas (Antarctic and Greenland) , which appear to generate 
artificial high pressure systems. This results in a "squeezing" of 
the pressure systems in lower latitudes, with consequent meridional 
structure quite different than the one in Figure 12, For T850 and 
G500 the annual range of 6,0 is overestimated by the model. 

From the annual cycles of the component 1,1 (Figures 24 and 25) 
it can be observed that the T850 cycle is overestimated in amplitude 
by the model, but with nearly correct phase. On the other hand, the 
model underestimates the amplitude of the 1,1 harmonic of G500, but 
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Figure 22. 
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Figure 24, . * 

~~ TIME SPECTRA FOR THE COMPONENT 'Ll 
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riqure 25, 

TIME SPECTRA FOR THE COMPONENT 
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rather faithfully reproduces the global phase. It also produces reason- 
ably close simulations of the 1,1 cycles of SLF. The Northern Hemisphere 
phase change (90° to 270°) from winter to summer (Figure 25) largely re- 
flects the seasonal change in pressure over Asia. 

The 2,1 harmonic (Figure 26) for the SLP represents mainly the dif- 
ference between the Siberian anticyclone and the Canadian high in the 

o 

Northern Hemisphere winter, as indicated by the 90 phase angle. The 
annual cycle of 2,1 is similar to 1,1, with an overestimation of the am- 
plitude, notably in July, and a seasonal phase reversal reflecting the 
monsoonal pressure change over Asia. In T850 the 2,1 largely reflects 
the seasonal temperature differences between the hemispheres as well as 
the heating of the Asiatic continent in summer, as indicated by the 
m aximum amplitude in July. The model correctly simulates the seasonal 
phase shift, but exaggerates the July amplitude, apparently due to over- 
heating of the Asiatic continent. For G500 the observed cycle of 2,1 
is generally followed in phase and amplitude by the model. 

In the SLP the cycle of 3,1 (Figures 27 and 28) is also dominated 
by the Siberian anticyclone during the winter, which is replaced by a 
low pressure system during summer. This is indicated, as in 2,1 and 
'1,1, by the seasonal phase shift. The above cycle is poorly represented 
by the model in the global analysis, while in the Northern Hemisphere 
the cycle is correct but overestimated, indicating some difficulty with 
the simulation of the Southern Hemisphere. The annual cycle of 3,1 for 
T85G is biased in amplitude for the Northern Hemisphere, while for the 
globe the summer cycle is not correctly represented in amplitude and 
phase . 

The component 2,2 (Figures 29 and 30) is related to the land- sea 
differences of temperature, geopotential and pressure. For all C fields 
similar annual cycles are observed, with two maxima and two minima. The 
model in July overestimates the amplitude of 2,2 in T850 and SLP, but 
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Figure 29. 


TIME SPECTRA FOR THE COMPONENT 
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the geographical positions of high and low values, as indicated by the 
phases, are well simulated. 

The 3,2 component (Figure 31) reflects not only the longitudinal 
land- sea differences, but also the asymmetry between the hemispheres. 

The same characteristics are observed in the annual cycle of 3,2 for 
T850 and SLP as in 2,2. The G500 cycle is poorly simulated by the 
model. This is due to the zonal structure of the above field, which 
is well accounted for by the zonal harmonics, leaving only a small 
perturbation to be represented by the higher harmonics. Thus 3,2 is a 
relatively minor harmonic for G500. 

The 4,2 component (Figures 32 and 33) ranks as a significant har- 
monic for all three fields. Additionally,' for G500 4,2 corresponds to 
the zonal distribution of the major ridges and troughs of the quasi- 
stationary long wave pattern. The Northern Hemisphere annual cycles 
of 4,2 are similar to the global cycles, but with amplitudes twice as 
large. This indicates that either the Northern Hemisphere is a good 
reflection of the Southern Hemisphere or, more likely, the characteris- 
tics of the Northern Hemisphere dominate for all fields. As in the 3,2 
harmonic, the M 4,2 component roughly parallels the annual cycle of the 
corresponding component for the actual climatology, with an overestima- 
tion of SLP and T850 in summer. The G500 cycle is fairly well simulated 
by the model in phases and amplitudes. 

Very little additional information about the model simulations can 
be derived from the annual cycles of other spectral components. This 
is apparent from the mismatch of the M and C components for n^4 and 
m§2. These harmonics representing smaller scale features of the syn- 
optic climatological patterns explain only a small part of the total 
variance of the three’ fields examined. However, it is possible to ob- 
tain further insight into the problem by comparing the observed and 
model climatologies as they are represented by their leading harmonics. 


- 33 - 



Figure 31 * 
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Figure 33» 


TIME SPECTRA FOR THE COMPONENT 
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This permits an objective comparison to be made between the large-scale 
features of the two climatologies through the use of correlation analy- 
sis. 

The following calculations were carried out: 

1. Seasonal mean fields of T850, G500, and SLP were reproduced ■ 
for C and M from their leading harmonics, ranked in descending order of 
amplitude, with the components accumulating and correlation coefficients 
were computed between the corresponding reproduced fields. 

2. A similar calculation was carried out using seasonal mean 
fields reproduced from harmonics ranked in ascending order of n and m. 

In all cases, the mean value 0,0 was omitted, and, in the case of 
T850 and G500, the meridional structure term 2,0 was also taken out be- 
. cause it results in an artificially large correlation coefficient be- 
tween the fields. The results obtained with both methods are displayed 
in Tables 9, 10, and 11 for T850, G500, and SLP respectively. 

From the tables it is observed that, as the harmonics used to. re- 
produce the fields accumulate, the correlation coefficient, r, either 
decreases or increases, depending on whether the fields were reproduced 
from the same or different dominant harmonics. For example, in Table 9, 
which lists the M and C T850 correlations, it is seen that when both 
fields are reproduced from the same single leading harmonic, the corre- 
lation coefficient is one. In spring when a second (mismatched) har- 
monic is added, the coefficient drops to 0.60, but when the third har- 
monics, which make up a set of three matching components, are added, 
the coefficient rises- up to 0.99. The addition of weaker harmonics 
then degrades the correlation. The summer and winter correlations of 
reproduced T850 fields are quite high, while the spring and fall values 
are much lower. In Table 10, however, it is seen that the correlations 
between the fields of G500 reproduced from the leading harmonics are 
quite low in winter and spring, but high in summer and fall. The SLP 
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Correlation coefficients between the reproduced seasonal M and C SLP 
fields from their; a) leading harmonics, b) harmonics ranked in 
ascending order of wavenumbers n,m. 
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correlations, shown in Table 11, are universally low, as expected, 
again indicating that the model fails to reproduce this field ade- 
quately. 
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Summary and Conclusions 


Spherical harmonic analysis has been applied to global and North- 
ern Hemispheric monthly mean observed (C) and model-generated (M) clima- 
tological fields of 850 mb temperature (T850) , 500 mb geopotential 
height (G500) , and sea- level pressure (SLP) for the purpose of evaluat- 
ing the atmospheric simulations produced by the GISS climate model. 

All three fields were found to be adequately represented by wave num- 
bers n,m^ 8, and, indeed, by a relatively small group of dominant har- 
monics for each field. Thus, the method of spherical harmonic analysis 

o 

does appear to be applicable, at least to the coarse-mesh (8 latitude 

o 

by 10 longitude) data sets employed. 

The C and M fields were first compared by ranking the dominant 
harmonics of both for each season in order of decreasing amplitude. 

The 0,0 harmonics, representing the global mean values, were found to 
be generally similar for the corresponding C and M fields. The model 
also successfully simulates the dominant 2,0 harmonics of T850 and G500, 
which represent the mean meridional gradients (i.e., the differences 
between Poles and Equator), although it does overestimate the amplitude 
of the temperature gradient as well as the global mean value, appar- 
ently due to overheating of the continents in summer. On the other 
hand, the model simulation of SLP is less satisfactory in terms of the 
leading harmonics. The 1,0 harmonics, representing the mean pressure 
difference between the hemispheres, are opposite in phase for C and M 
summer, indicating an incorrect mass distribution, while the 2,0 har- 
monics reflect a tendency of the model to build high pressure over the 
Antarctic and low pressure at the Equator. 

Monthly spherical harmonic analyses for both the globe and the 
Northern Hemisphere were obtained for the M and C fields. Then the 
annual cycles of the C leading components were compared with the 
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corresponding M harmonics in terms of amplitude and phase differences. 

The annual cycles of the global mean values were generally correctly 

simulated by the model, while for the Northern Hemisphere the hemispheric 

mean values for T850 were overestimated in summer. Also, the hemispheric 

means for SLP were low biased for all months. This was reflected in the 

o 

M annual cycle of the 1,0 harmonic, with a phase shift of 180 from 
April through August, indicating, as before, an incorrect pressure dis- 
tribution between the hemispheres. The M 2,0 annual cycles of T850 
and G500 were parallel to those of C, while for SLP the same term was 
incorrectly simulated in phase and amplitude by the model. The M 4,0 
annual cycles of T850 and G500 were overestimated in simmer, especially 
for the Northern Hemisphere, due to the overheating of the continents. 
This becomes more apparent from the annual cycles of the M 2,2, 3,2, 4,2 
terms, which represent the continent -ocean temperature differences for 
T850 (overestimated in summer) . The simulated temperature distributions 
resulted in overestimations (underestimations) in the amplitude of the 
zonal pressure distributions for summer (winter) . This was interpreted 
as a failure of the model to simulate the four basic pressure systems 
of the Northern Hemisphere with the correct intensities. The annual 
cycle of the M meridional pressure distribution, 6,0, was not similar 
to the climatological one (Figure 12) . This indicates the false meri- 
dional structure of the predicted SLP field. 

An objective comparison of the large-scale features of the two 
climatologies was obtained by correlating the reproduced seasonal M and 
C fields from their components * The results taken from this method 
generally confirm the subjective comparison of the C and M leading har- 
monics for the same fields. Also, they indicate the fact that G500 and 
T850 were successfully simulated by the model (especially for the sum- 
mer) , while the SLP was poorly simulated for all seasons . 
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. • APPENDIX A. 


Normalization of the spherical harmonic functions 


If P (sin<j>) are the unnormalized Legendre associated functions, and 

XI fill 

- jn P (simfO 

P n,m (s:Ln<i,) = (cos< ^ *d{sin<j>) m ' 

where P n (sin<}>) are the Legendre polynomials defined as 
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Thus the complete normalizing factor for this function, in order for it's mean 
square value over a sphere to be unity, is 

[ 2(2n-KL) (n-m) 1 ~| ^ 

S m (n4m) ! J 

It) normalize the sperical harmonic function to (4m)" 1 , we use the following 
set of harmonic functions: 
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APPENDIX B 


Evaluation of the expansion coefficients . 


Die observed field, Q(<j>,X), is represented by the series. 


oV-x) 


= 11 

n,m 


{ C *Y (el + S ) 

n,m n,m n,m n,m 


Multiplying both sides by ^ and integrating over the domain of 0(<}>,A) , (i.e. 
the -domain of the observations of Q(<j>,X) as gridded) , we get. 




( C *Y fe ^ »Y^\ 
n,m n,xn i,j-. 


+ S *Y W *Y fk \ 
n,m n,m 1,3 


)*dft 


where dfl — cos<f>d({)dX and tt/2 < $< tt/2 , 0 ^ X^ 2 ir and s is a spherical surface. 

If the superscript (JO indicates evenness of the Y , it follows that 

n,m 



1 1 

n,m 
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, yfe^ 

n,m i, j 


dfi 


+ 0. 


Expanding the sumnations and taking into account the orthonormal ity property of the 
spherical harmonic functions. 
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C. . *1 and 
i#3 



J 3 
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An analogous relation is derived for S. . 
Finally, 


if we let (k'l denote oddness of the Y 

n,m 


n,m 


n,m. 


27J> % 


■Ui 



QOM) 


Y (e1 , 
n,m 

Jo) 

n,irt 


•cos<|> dtp dX 


If Q(4>,A) = Q(~cf>,A) , which states that Q is an even function of latitude <p. 


then 


( C n,rai ,. 27T cosmX to; 

J a. * •/« 

n,m J 0 sinmA /_t fe 


TVZ 

Q(<f>,A) »P n m (sin<j)) d4> 


vhere N is the normalization factor. Uie second integral above has symmetrical limits, 
for (n-m) — odd , the function P (sin<f>) is odd, cos<j> is even and therefore the 

11 /Iu 

quantity Q(<p,X) *P n ^(sin^) *cos<}) is an odd function of latitude <p, and 



n,m 
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If the data Q are also functions of time (form a time series) , then 


[W w i 


S n, m <« 


= j QM,A,t)-j 
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When averaged over a tine period T, the last equation becomes. 
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Y fe) 
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where ( ) 


/ 3( > dt 

T i 


and T. - T . = T 

] X 


2hus we may either average the fields Q over time, and apply spherical harmonic analysi 

to Q, or apply the harmonic analysis to each field Q and average the expansion 

coefficients/ r . C 

n,m ,n r m 
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